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T-H I I discuss how masses and widths of hadron resonances are extracted 

^ . from lattice QCD. Recent lattice results on the light, strange and charm 

CO ' meson resonances are reviewed. Their properties are revealed by simulat- 

'^ , ing the corresponding scattering channels tttt, Ktt and Dn on the lattice 

^ ■ and extracting the scattering phase shifts. In particular we address the 

^ ; resonances p, i:';5(2400), i:'i(2430), K*, k and i^^(1430). 
■^ . 

en ' 1. Introduction 



X 



Most hadrons are resonances, i.e, they decay extremely fast via the 
strong interaction. Yet most of these resonances were studied in lattice 
QCD assuming the so-called narrow width approximation, that is ignoring 
C^ , their strong decay. Up to now only the p meson has been simulated properly 

as a resonance by several groups and its width was extracted (see [H [21 [3] 
and references therein) . Recently the first simulation of charmed resonances 
in Dtt scattering was performed [1] , while the strange scalar and vector res- 
onances were addressed by simulating Ktt scattering [S]. This paper briefly 
reviews the main results and methods that were employed. Recent lattice 
results of resonances are reviewed in [6] . 
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2. Meson-meson scattering in a resonant channel on the lattice 

Meson resonances are formed in the strong scattering M1M2 — )• i? — )• 
M1M2 in partial wave /. They exhibit a Breit-Wigner-hke resonance behav- 
ior of 6, amplitude T and a{s) oc sin^ 6{s) 



(p*)2'+i cot 5 



2 v^ - ml) (2) 



vs g 

where r(s) is parametrized in terms of the phase space and the R — )■ M1M2 
coupling g. The combination (p*) '^^ cot 6 / ^/s is linear in s for a single 
Breit-Wigner resonance, which allows the extraction of mR and g (and 
therefore the width) using a linear fit ([2]) once the phase shifts 6{s) are 
determined from the lattice. So the goal is to simulate the scattering on the 
lattice and determine the scattering phase shift S{s). 

For this purpose one computes the correlator Cij{t) = (0104(^)0^(0)10). 
We use the interpolators O = Mi{pi)M2{p2) = Qi^iQi Q2^2Q2 that create 
two-meson states with definite momenta, and O = qTq' that couple well 
to the resonances. Both are constructed to have the quantum numberq^ of 
the desired channel and total momentum P = pi+P2- The interpolators O 
couple in general to all physical eigenstates n and each of them evolves as 
^-E„t j^ euclidean time, so Cij{t) = ^nAf e"-^"*. 

We calculate the correlators Cij{t) using a powerful distillation method 
[7], which enables the calculation of all the necessary Wick contractions. 
Our study is based on 280 gauge configurations with a ~ 0.124 fm and 
dynamical Wilson-Clover u/d quarks corresponding to m-^ ~ 266 MeV. A 
rather small volume 16^ x 32 makes the distillation method [7] feasible. The 
valence charm quark is treated using the Fermilab method described in [Ij. 
The resulting Cij{t) allows the extraction of the few lowest eigen energies 
En via the generalized eigenvalue method. 

3. Physics information based on the energy spectrum En{L) 

The scattering energy levels (black and green in Fig. [1]) appear at 
E{L) = \/m\ -I- p^ -I- \/m^ -|- p^ + l^E{L) with discrete pi = n^ due to the 
periodic boundary conditions in space. The energy shift AE(L) in the finite 
volume is due to the strong interaction of the two mesons. The negative 
shift of the lowest level for / = 1/2 s-wave scattering of Kit, Dtt and D*Tr 



^ On a finite discrete lattice the interpolators O have to transform according to irre- 
ducible representations of the symmetry group related to the center-of-momentum 
frame. 
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in Fig. [1] indicates attractive interaction. The positive shift for Kir with 
/ = 3/2 indicates repulsive interaction. In addition the the scattering levels 
shown in black and green, the presence of resonances in I = 1/2 channels 
leads to the levels drawn in red and violet. These indicate s-wave resonances 
L>o(2400) in Dvr, L>i(2420) and L>i(2430) in L>*7r, Kq{U^Q) in Kvr, as well 
as p-wave resonances K*{^92), K*(1410), K*(1680) in /Cvr [11[5]. We do 
not find an additional level due to k [5], which is in agreement with the 
fact that the experimental phase shift does not reach 90° below 1 GeV. We 
also do not find additional energy levels in I = 3/2 channels, in line with 
absence of the exotic resonances in experiment. 

Kti, s-wave, 1=1/2 s-wave, 1=3/2 Kti, p-wave, 1=1/2 p-wave, 1=3/2 
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Fig. 1. The energy levels (effective masses of eigenvalues) for Kit [J^ = 0+, 1 ), 
D-K (J^ = 0+) and D*it (J^ = 1+) scattering with P = i]. Dashed lines 
indicate energies of non-interacting scattering states. 



4. Phase shifts and resonance parameters 

The energy shift in finite volume reveals the attractive or repulsive na- 
ture of the interaction. However, it also rigorously renders the phase shift 
for the elastic scattering in the infinite volume via the Liischer's relation. 
In particular, the energy level E{L) for a scattering system with momenta 
P renders the elastic phase shift 5{s) at s = E'^ — P'^ in partial wave /. Note 
that the extraction of 5{s) is straightforward only when the partial-wave 
mixing due to the discrete symmetry is absent or negligible, which usually 
holds well for P = 0, but holds rarely for the scattering of two particles 
with different mass and P ^ [S]. 
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Fig. 2. Left: the p-wave tttt — ^ tttt phase shift 6 with / = 1 [1 . Right: the 
corresponding {p*)^ cotS/^/s ^. 

The TTTT — )• p — )• TTTT was simulated for three different P in [1] and five 
energy levels lead to resonant phase shift in Fig. [21 The linear fit of the 



resulting p*^cot5/^/s <^ leads to g^^^^ = VOttq = 5.13 ± 0.20 and m^"* = 
792 ± 10 MeV compared to gp^n = 5.97 and rrip^'P = 775 MeV. This is the 
only resonance where proper lattice treatment has reached a certain level of 
maturity (see [H El [3] and references therein). A particularly detailed and 
impressive shape of the resonant phase shift curve was achieved in [3j at a 
heavier pion mass mj^ ~ 400 MeV. 

Over the past year we performed the first simulation of the Kit ^, Dtt, 
D*TT [1] and pir [9] scattering and the corresponding resonances. Since this 
involves scattering of two particles with different masses, we considered only 
P = when the mixing of different I is absent or negligible |Sj • 

The energy levels for Kir scattering in Fig. [Ulead to the phase shifts in 
Fig. [31 for s-wave and p-wave with I = 1/2, 3/2. These are in qualitative 
agreement with the experimental ones. 

Next we concentrate on the charmed resonances that appear in Dtt and 
D*7r. The three energy levels for Dtt s-wave scattering in Fig. [H lead to 
the phase shifts [3], and the linear fit ^ over three points leads to m and 
r (or g) for the broad Z)g(2400) in Table [H The analysis of D*7r spectrum 
with J^ = l"*" is more complicated since there are two nearby resonances in 
experiment, as evidenced also by the red and violet levels in Fig. [H We find 
that the red level for D*Tr scattering is due to the narrow Di(2420) which 
decays only in d-wave in the rUc — t- oo limit. In this limit the remaining 
three levels are related to s-wave D*7r scattering which is dominated by the 
broad Z)i(2430). A linear fit ([2]) through these three points leads to the 
mass and the width for the broad L'i(2430) in Table [H 

The resulting masses and widths of Dq (2400) and Di (2430) agree quite 
well with the experimental ones. Since Dq(2400), located at ~ 2318 MeV, 
is very close to its strange partner 15*0(2317), several authors proposed 
that Dq{2A00) has a sizable tetraquark component cssu. We get w.£)*(2400) 
near the experimental value without explicitly incorporating the additional 
strange valence paio 



The ss can not appear as intermediate state in our simulation without dynamical s. 
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Fig. 3. The Ktt phase shifts 5l in channels / = 0, 1 and / — 1/2, 3/2 as a function 
of Kit invariant mass y/s. The lattice results are given by the red circles (they 
apply for ttItt — 266 MeV) [5], while the other points are experimental phase shifts. 



"T-D^(2400) - Trn g D*(2Am)^DTT 



"^Di(2430) ~ ""^ 9 Dx{2im)^D-K 



lat[l] 
exp 



351 ± 21 MeV 2.55 ± 0.21 GeV 
347 ± 29 MeV 1.92 ± 0.14 GeV 



381 ± 20 MeV 2.01 ± 0.15 GeV 
456 ± 40 MeV 2.50 ± 40 GeV 



Table 1. The charmed resonance masses (with respect to m = ™-D+^™g* ^ ^j^fj ^]^g 



VjV^U UV^ Iff, . 

couplings 5, which parametrize the widths F = g^p* / s. The experimental couplings 
g are derived from total widths. 



The compilation of the D meson spectrum in Fig. S] shows quite good 
agreement with experiment [1]. The masses of broad resonances Dq(2400) 
and Di(2430) are extracted as explained above. Other four low-lying J^ = 
0^, 1", 1"*", 2+ states are stable or very narrow, so they were simulated using 
O = cTu and ra = E (P = 0) is employed like in all previous simulations. 
This narrow-width approximation is applied also for the excited states in 
J = 0~, 1~, 2~ channels, which are compared to the states observed by 
BaBar in 2010 [10] ; unfortunately these are not yet confirmed by any other 
experiment. 

The scattering lengths for i^vr, Dtt and 0*71 were also extracted in [5l[3]. 

In conclusion, p is the only resonance that was treated properly by sev- 
eral lattice groups up to now. We presented the first results of the strange 
and charmed resonances based on the simulation of the corresponding scat- 
tering channels. 
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Fig. 4. Energy differences m — |(m.D + Sto^.) for 13 mesons on lattice [1] and in 
experiment; the reference mass is to = j{mD + Sm^)*) w 1971 MeV in experiment. 
Magenta diamonds give masses for states simulated as resonances [1]. Masses 
extracted as energy levels on a finite lattice are displayed as blue crosses [i]. 
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